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^ , Abstract 

^ ' We apply a recently suggested technique of the Neumann-Dirichlet reduction to a toy model 

, of brane-induced gravity for the calculation of its quantum one-loop effective action. This model 

is represented by a massive scalar field in the (d -I- l)-dimensional fiat bulk supplied with the 
d-dimensional kinetic term localized on a fiat brane and mimicking the brane Einstein term of the 
Dvali-Gabadadze-Porrati (DGP) model. We obtain the inverse mass expansion of the effective 
action and its ultraviolet divergences which turn out to be non-vanishing for both even and odd 
spacetime dimensionality d. For the massless case, which corresponds to a limit of the toy DGP 
■ model, we obtain the Coleman- Weinberg type effective potential of the system. We also obtain the 

' proper time expansion of the heat kernel in this model associated with the generalized Neumann 

Qf^ , boundary conditions containing second order tangential derivatives. We show that in addition to 

T-H ' the usual integer and half-integer powers of the proper time this expansion exhibits, depending 

T-H ' on the dimension d, either logarithmic terms or powers multiple of one quarter. This property 

, is considered in the context of strong ellipticity of the boundary value problem, which can be 

^— «' ' violated when the Euclidean action of the theory is not positive definite. 

^: 

p • 1. Introduction 

(D : 

I Modified theories of gravity in the form of braneworld models can in principle account for the phe- 

^ I nomenon of dark energy as well as for nontrivial compactifications of multi-dimensional string models. 

It becomes increasingly more obvious that one should include in such models the analysis of quantum 
effects beyond the tree- level approximation [T]. This is the only way to reach an ultimate conclu- 
. sion on the resolution of such problems as the presence of ghosts [2] and low strong-coupling scale 

" " ' [3] . Quantum effects in brane models are also important for the stabilization of extra dimensions [1] , 

fixing the cross-over scale in the Brans-Dicke modification of the DGP model |5j and in the recently 
suggested mechanism of the cosmological acceleration generated by the four-dimensional conformal 
anomaly [S]. 

A general framework for treating quantum effective actions in brane models (or, more generally, 
models with timelike and spacelike boundaries) was recently suggested in [71 [HI [5] ■ The main pecu- 
liarity of these models is that due to quantum field fiuctuations on the branes the field propagator 
is subject to generalized Neumann boundary conditions involving normal and tangential derivatives 
on the brane/boundary surfaces. This presents both technical and conceptual difficulties, because 
such boundary conditions are much harder to handle than the simple Dirichlet ones. The method 
of [9] provides a systematic reduction of the generalized Neumann boundary conditions to Dirichlet 
conditions. As a byproduct it disentangles from the quantum effective action the contribution of the 
surface modes mediating the brane-to-brane propagation (that is, within one brane), which play a 
very important role in the zero- mode localization mechanism of the Randall-Sundrum type [TU] • The 
main purpose of our paper here is to apply this method to a simplified toy model of brane-induced 
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gravity in order to demonstrate how it works for the first nontrivial example of a field system with 
boundary conditions involving second-order tangential derivatives. As we will see, this model leads to 
qualitatively new structures in the one-loop effective action, its renormalization, and the associated 
heat kernel. 

Briefly the method of [9] looks as follows. The action of a (free field) brane model generally contains 
the bulk and the brane parts, 

SiH^l; I d^+^Xcj,{X)F{V)<^{X) + \ [ d''x^ix)^{d)ip{x) , (1.1) 



2 Jb 2 jj, 

where the (d + l)-dimensional bulk and the d-dimensional branc coordinates are labeled respectively 
hy X = X^ and x = x'^, and the boundary values of bulk fields (j>{X) on the brane/boundary b — 9B 
are denoted by ip{x), 

0(X) ^^(x). (1.2) 

b 

The kernel of the bulk Lagrangian is given by the second order differential operator i^(V), whose 
derivatives V = dx are integrated by parts in such a way that they form bilinear combinations of first 
order derivatives acting on two different fields (this is denoted in (|l.ip by the double-headed arrow). 
Integration by parts in the bulk gives nontrivial surface terms on the brane/boundary. In particular, 
this operation results in the Wronskian relation for generic test functions 2(^)1 

j d''+^x(^cPiF{V)q^2-(biF{y)h)=-J^^d''x(^cj)iWh-(biW(t>2). (1.3) 

Arrows everywhere here indicate the direction of action of derivatives either on (jji or 02 • 

The brane part of the action contains as a kernel some local operator k(9), d — dx- Its order 
in derivatives depends on the model in question. In the Randall-Sundrum model [lOj . for example, 
it is for certain gauges just an ultralocal multiplication operator generated by the tension term on 
the brane. In the Dvali-Gabadadze-Porrati (DGP) model [11] this is a second order operator induced 
by the brane Einstein term on the brane, k(9) ~ 99/to, where m is the DGP scale which is of the 
order of magnitude of the horizon scale, being responsible for the cosmological acceleration [12]. In the 
context of the Born-Infeld action in D-brane string theory with vector gauge fields, K,{d) is a first-order 
operator [13]. 

In all these cases the variational procedure for the action (jl.ip with dynamical (not fixed) fields 
on the boundary ^{x) naturally leads to generalized Neumann boundary conditions of the form 

(vK(V) + ^(9))</)|^-0, (1.4) 

which uniquely specify the propagator of quantum fields and, therefore, a complete Feynman diagram- 
matic technique for the system in question. The method of [3] allows one to systematically reduce 
this diagrammatic technique to the one subject to the Dirichlet boundary conditions 0| b =0. The 
main additional ingredient of this reduction procedure is the brane operator F^^^^^''{x,x') which is 
constructed from the Dirichlet Green's function G£i{X^X') of the operator FiV) in the bulk, 

F'^™^^(x,x') = -M^(Vx)Gz5(X,X')V^(Vx') + K{d) 5{x,x') . (1.5) 

X—e{x),X'—e{x') 

This expression expresses the fact that the kernel of the Dirichlet Green's function is being acted 
upon both arguments by the Wronskian operators with a subsequent restriction to the brane, with 
X = e{x) denoting the brane embedding function. 

As shown in [3] , this operator determines the brane-to-brane propagation of the physical modes in 
the system with the classical action (|l.ip (its inverse is the brane-to-brane propagator) and additively 
contributes to its full one-loop effective action according to 

ri_ioop = \ Tr(v'+'^ ^'^P=\ + ^ InFi^-'-, (1.6) 
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where Tr^"'^' denotes functional traces of the bulk theory subject to Dirichlet and Neumann boundary 

conditions, respectively, while Tr^''' is a functional trace in the boundary d-dimensional theory. The 
full quantum effective action of this model is obviously given by the functional determinant of the 
operator i^(V) subject to the generalized Neumann boundary conditions (|1.5[) . and the above equation 
reduces its calculation to that of the Dirichlet boundary conditions plus the contribution of the brane- 
to-brane propagation. 

Here we apply (jl.6p to a simple model of a scalar field which mimicks in particular the proper- 
ties of the brane-induced gravity models and the DGP model TT]. This is the {d + l)-dimensional 
massive scalar field (j){X) = 4>{x, y) with mass AI living in the half-space y > Q with the additional 
c?-dimensional kinetic term for (f{x) = 4i{x, 0) localized at the brane (boundary) at y = 0, 

SiH^l J d''+^x({Vc}>{X)f + M^4>^{X))+^jd''x{dv{x)f. (1.7) 

Here and in what follows we work in a flat Euclidean (positive-signature) spacetime. Therefore, this 
action corresponds to the following choice of F{V) in terms of (rf+ l)-dimensional and d-dimensional 
D'Alembertians (Laplacians) 

F(V) = m2 - = m2 -U-dl, □ = = dl- (1.8) 

Its Wronskian operator is given by the normal derivative with respect to outward-pointing normals 
to the brane, W = —dy, and the boundary operator k(i9) equals 

Km) = -—a, (1.9) 

m 

where the dimensional parameter m mimicks the role of the DGP scale [11 . Thus, the generalized 
Neumann boundary conditions in this model involve second-order derivatives tangential to the brane, 



dy + -D 
m 



0, (1.10) 



cf. (fL4| with W = -dy and k = -D/m. 

As we show below, the brane-to-brane operator for such a model has the form of the following 
pseudodifferential operator, 

F^''''"'=(9) = — (-D + rnVM^^D). (1.11) 
m 

In the massless case of the DGP model [H], M = 0, this operator is known to mediate the gravitational 
interaction on the brane, interpolating between the four-dimcnsional Newtonian law at intermediate 
distances and the five-dimensional law at the horizon scale ~ 1/m [3]. We calculate the effective 
action p.6p for this model both in the form of the l/A//-expansion and exactly in terms of a special 
hypergeometric function representation (as a function of M and the dimensionless ratio m/M) and 
find its ultraviolet divergences. 

As a byproduct of the Neumann to Dirichlet reduction (|1.6p . the technique of ^9] also yields a 
method for obtaining the proper time expansion for the heat kernel associated with the boundary 
conditions (jl.4p . For simple Neumann (Robin) boundary conditions containing at most first order 
derivatives tangential to the boundary, this expansion has the form 



oo 

Tr^'^'^ e-^-^-^ = j^^^;^ E {^^ A,. + s"/^ B.„,,) . (1.12) 



n=Q 



In addition to the well-known bulk integrals An of the local Schwinger-DeWitt coefficients of inte- 
ger powers of the proper time [HI [TH [16], it contains surface integrals Bn/2 as coefficients of both 
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integer and half-integer powers of s [T71 [THl [111 [201 HI] ■ They are sufficiently easy to calculate for 
the Dirichlet boundary conditions [T71 [H], but become much harder to obtain for the Robin and 
generalized Neumann case with a growing number of tangential derivatives [THl HH US] • As shown in 
[9] , the Neumann-Dirichlet reduction method simplifies their calculation essentially. For second-order 
derivatives they are not known at all, and a toy DGP model of the above type seems to be the first 
application of the heat kernel mehtod subject to the boundary conditions (jl.lOp . 

It turns out that in the case of (jl.lO|) even the very structure ()1.12p is incorrect, because for even 
d it contains also terms logarithmic in s and for odd d it has also powers of s which multiples of a 
quarter. We calculate these additional terms in the heat kernel expansion, discuss their relation to 
a nontrivial analytic structure of the brane part of the effective action (|1.6p and also to the problem 
of strong ellipticity [23] of the boundary value problem (jl.lOp . In particular, the latter is shown to 
be determined by the positivity of the action (|1.7p or the positive-definiteness of the brane operator 

(iniD. 

The paper is organized as follows. In Sect. 2 we derive the Dirichlet and brane parts of the 
effective action (|1.6p . Sects. 3 and 4 present its inverse mass expansion and ultraviolet divergences in 
various dimensions. In Sect. 5 we obtain a new type of the proper time expansion for the heat kernel 
associated with the boundary conditions involving second-order tangential derivatives. In Sect. 6 we 
present the hypergeometric function representation of the effective action and analyze the limit of a 
simple Neumann boundary condition corresponding to to — > oo; we also consider the massless limit 
M = which gives the effective potential in the toy DGP model. In the Conclusions we discuss these 
results in the context of a possible violation of strong ellipticity for the boundary conditions (|1.10p . 
their potential applications in braneworld models including gravitation and the use of the proper time 
method in brane models. Three appendices contain the derivation of the inverse mass expansion of 
the effective action, its exact hypergeometric function representations, and the presentation of the 
status of the strong ellipticity problem in a toy DGP model. 



2. Dirichlet and brane-to-brane contributions 

We begin by constructing the Dirichlet part of the effective action. The basic Dirichlet Green's function 
of the model can be obtained by the proper time integration of the corresponding heat equation kernel, 

Kd{s\X,X') = e'°(''+^i5{X,X'). (2.1) 

It follows by the method of images from its well-known expression in the flat {d + l)-dimcnsional 
spacetime without boundaries 

Kd(s |A,X ) = — — -j^ \ exp 



(47rs)" 



4s 

-exp(-^^'-')'+/^ + ^')") ). (2.2) 



The functional trace of this heat kernel contains two terms — bulk and boundary integrals of the 
only two nonvanishing Schwinger-DeWitt coefficients, 



( d'^+^X Kd{s\X,X) 
Jy>a 



'([ d^^^X-^s^'- fd'^x). 

S 2 \Jy>Q ^ J J 



(2.3) 



(47rs) 2 \Jy>0 

The corresponding Dirichlet-type effective action for the model with the mass M in the bulk can be 
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obtained by the following proper time integration, 



i Tr_D In 



(d+i) - 
An 



1 



ds 



2 . n S 



r - 



-s Af 

47r 



(2.4) 



The dimensionality d will be treated as a parameter of the dimensional regularization. Therefore, this 
expression contains ultraviolet divergences as the poles of Gamma functions at negative integer values 
of their arguments. These divergences are represented here either by the bulk or boundary surface 
integrals, depending on whether the total spacetime dimensionality {d + 1) is even or odd. 

For the brane part of the action (|1.6p we need the brane-to-brane operator (|1.5p which is based 
on the Dirichlct Green's function of the model. The latter is also exactly calculable in elementary 
functions because the corresponding proper time integral can be expressed in terms of the modified 
Bessel function of half-integer order. 



Gd{X,X')^ J dsKD{s\X,X 



ds 
1 



-\y-y'\ VA/^-D 



-{y+y')VM'^-n 



S{x, x') 



- □ 

Therefore the first term of (jl.Sp takes the form of a square-root operator [5] 



dyGD{X,X')dy 



X={xfi),Y=[x' fi) 



= v/M2 -U5{x,x') 



(2.5) 



(2.6) 



and the full brane-to-brane operator (jl.Sp is given by (|l.lip . 

The operator (|l.lip is of a nonlocal pseudodifferential nature, and no conventional proper time 
representation is known for its functional determinant (see Sect. 7, though). Therefore we will calculate 



the latter in the basis of Fourier modes — the eigenmodes of F 
delta-function in the Fourier integral we have 

- Tr^-^MnF^™'" ^ - f d' 



branc 



By resolving the d-dimensional 



In (^-a + m\/ - S{x,x') 



d'^a 



{inrPrid/2) J JO 
where p is the radial integration variable in the momentum space 



dp/"^ In + my^Ap+p^ ) 



(2.7) 



(2. 



As we see, the mass parameter enters the logarithmic function here in a very nontrivial way, so 
a typical l/M^-expansion of the local Schwinger-DeWitt expansion is far from being straightforward. 
In the next section we derive this expansion by converting the expression (j2.7p into the form of a 
so-called integral with a weak singularity to which a known asymptotic expansion technique can be 
directly applied. 
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3. Inverse mass expansion 



By integrating in ()2.7|) by parts and using the rules of the dimensional regularization, which annihilates 
purely power-divergent integrals, we have 



J dpp'^-^ In [jp- + to\/m2 





_^2 + m/^W+ 







pi 



Then, with the change of the integration variable from p to 

„2 r -,1/2 



t = 



P 



P : 



M 



p^ 



2t{Vl + i^ + t) 



the integral takes the form 

dpp'^^^ In (^p^ + m^jNP + p'^ 



I- 





2 - 



(2M2)''/2 



where 



/ = 



dtf^'^-^it + ey^^it), 



(Vl + t2+i)'i/2, 



and 



m 

2M' 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
(3.5) 

(3.6) 



Obviously the 1/M asymptotic expansion corresponds to the asymptotic expansion in e ^ — 
the limit in which the power-law singularity of the integrand occurs at the lower integration limit 
t = 0. Here the integrand is not analytic because of the factor t'^/^-i having a branch point at 
t = 0. Remember that the integral should be calculated for a generic dimensionality d that should be 
analytically continued to the complex plane in order to regularize the ultraviolet divergences appearing 
at the upper integration limit t — s- oo. Therefore, one should expect that the expansion of the integral 
will also have a part nonanalytic in e — > 0. 

As shown in Appendix A by the asymptotic expansion method for integrals with a weak singularity 
[25] ■ this expansion indeed has the form 



d/2-l 



i(7rd/2) 



.2j 



(3.7) 



where the first term — the nonanalytic part — has a branch-point singularity. Interestingly, the 
coefficient of this nonanalytic factor is exactly expressed through the function (|3.5p itself but with the 
flipped sign of the argument. Its expansion in higher powers of e, given by Eq. (|A.9p in Appendix A, 
is determined by the derivatives (^'■"-'(0) = d^'^ (p / dV^ (0) which explicitly equal 



^(")(0) 



d fd 



r(i 



2 V 2 ^ 



2 - (" - 4) 



The analytic part of (|3.7p contains only even powers of e, with the coefficients 

1 - 2j) 



a2j 



^ _22j-^-d/2 ^ r(i V-)r(| 



(3.8) 



(3.9) 
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Thus, finally, on account of (|2.7p and (|3.3p the inverse mass expansion for the brane part of 
the action takes the form of two terms having qualitatively different analytic behaviour in the mass 
parameters M and m, 



- Tr*-"^' In F '^'^'^"^ 
2 



dT{d/2) 



f — I f{-e) + 



An 



27r 



(3.10) 



4. Ultraviolet divergences versus spurious infrared poles 



The integral (|3.4p has a power divergence of order d—1 at the upper integration limit. Its differentiation 
with respect to e improves the convergence of the integral which becomes ultraviolet finite after d 
differentiations. This means that the ultraviolet divergence of p.4p is a polynomial in e of order d—1. 
The ultraviolet divergence of the brane action is a polynomial of order d in m and M, respectively, 
cf. Eq. p.3p . In dimensional regularization these divergences manifest themselves as poles in the 
dimensionality d analytically continued to the complex plane. From (|3.10p it follows, however, that 
for even d the nonanalytic and analytic parts of the inverse mass expansion separately have poles 
to all orders in e — the poles of tt/ sm{7Td/2) and the poles of one of the Gamma functions in the 
numerator of p.9p for all 2j > d/2 — 1. This implies an intrinsic cancelation between the infinite 
sequence of poles in the nonanalytic and analytic parts of (|3.10p . This cancelation of spurious poles, 
which have the nature of infrared divergences, can be directly observed by calculating separately these 
two contributions. 

The real ultraviolet divergences of (j3.4p can be independently obtained by means of integration by 
parts. In terms of the new integration variable x = 1/t the integral becomes divergent at the lower 
integration limit, 

1= / dxx-''f{x), (4.1) 
Jo 

fix) = (l + £a:)-i(Vl + .x2 + l)l (4.2) 

We call in the following the physical dimension of spacetime N in order to distinguish it from the 
formal dimension used in the integrals. Analytically continuing the spacetime dimensionality from 
its physical value N to the domain of convergence, where integrations by parts are possible without 
introducing extra surface terms, we have 

Then taking the limit to the real physical dimensionality, 

d = N -S, ^ ^ +0, (4.4) 
we have the logarithmic divergence of this integral as a residue of the pole in 5 0, 

S r(N) ^ ' 

By applying this formula we obtain the ultraviolet divergences which for even and odd dimension- 
alities look as follows. For even d being a multiple of 4, d = 4fc — 5, fc = 0, 1, • • • , the divergences 
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read 



- Tr^''^ In F 



5 2{AnYnT (I) ] ""^^^ r(2fc - j)(2j + 2 - 2k)\ ' ^^'^^ 
For even d being a multiple of 2, c? = 4fc + 2 — i5, fc = 0, 1, • • • , they have the form 



1 



Tr^'^) InF 



raiie 



div 



^ 2(47r)'i/2r (f ) 7 ''^^r(2fc + l-j)(2j + l-2fc)!- ^''^ 

Similarly, for odd d = 2fc + 1 — (5 the ultraviolet divergences of the brane effective action read 

1 ( j\ I 'ii'^ 

- Tr*- In F 
2 

= i X i / d^'^+lx V ("l)^^^^^r(fc-2j- 2fc-2j 2J + 1 .4 

^"2(4.)'^/2r(f)r r(i-j)(fc-j)! • ^'-'^ 

These results can be directly checked by disentangling the poles in spacetime dimensionality in 
the general inverse mass expansion p.lOp . This calculation confirms the cancelation of the fictitious 
infrared divergences in even d mentioned above. Although these divergences have a spurious nature and 
only arise at intermediate calculational stages, their presence reflects a nontrivial analytic structure of 
the asymptotic mass expansion p.lOp . As we will see below, they entail a nontrivial form of the heat 
kernel expansion corresponding to the generalized Neumann boundary conditions with second-order 
derivatives tangential to the boundary. 



5. Heat kernel expansion 



Because of the nontrivial pseudodifferential nature of the brane-to-brane operator (jl.lip and the 
inverse mass expansion of its action (|3.10p the heat kernel for the generalized Neumann boundary 
conditions (jl.lOp docs not have a typical expansion in integer and half-integer powers of the proper 
time ()1.12p . Thus we assume a more general structure of this expansion in the form 



(47rs) 



ly>0 



{P} 



d'^xb. 



{p} 



(5.1) 



where the summation runs over some unknown set of powers {p} of the proper time with some unknown 
coefficients of the surface integrals b^py . The bulk integral here involves only one term corresponding 
to the only non- vanishing bulk Schwinger-DeWitt coefficient aQ{X,X) = 1 (which is independent of 
the boundary conditions and coincides with the one in spacetime without boundaries). 

This expansion generates the inverse mass expansion for the effective action of the problem (jl.ip . 



i Trjv In 



- □ 



(d+i) 



^ Tin I — e 



I y>0 



{p} 



M2P 



d'^xb 



{p} 



(5.2) 
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Our goal now will be to determine the range of summation {p} and the coefficients 6{p} by comparing 
this expression with the inverse mass expansion of (|1.6p known from (|2.4p and (|3.10p . 

We begin by considering the most interesting case of even dimensionality which under dimensional 
regularization reads as d = 2fc — ^, ^ — 0. Even though the heat kernel is an ultraviolet finite object 
(which generates UV divergences in the action due to the divergence of the proper time integration) , 
we need this regularization to regulate intermediate infrared divergences which cancel out in the final 
answer. Thus, assembling together p.4p and (|3.10|) we get 

r 



1 



■ Tr jv In 



- □ 



(d+i) 



An 



1 



-id+l 



X 



y>0 



1 



47r 



r(-d/2) 



E 

n—k 



-r(-d/2) + X^^;^a2, 



,.0 ^r(f) 

1 (^("-'^•)(0) to"-*/2 



{n-k)\ M"-*V2 



d^2k-S. 



(5.3) 



Comparing it with (|5.2|) we immediately get the range of summation over the proper-time powers {p} 
which we will label by the integer numbers j and n, 



ri/2, 

M = < J, J > 1, 

[n/2-5/A, n>k + l, 
and the corresponding coefficients 



(5.4) 



Hp} 




J > 1, 



bn/2-S/i = K/2, n> k+1. 
In terms of 61/21 bj and 6„/2 the heat kernel trace takes the form 



Tr 



1 



N 



(Atts) 



y>a 



(5.5) 



(47rs)" 



61/2 s 



1/2 



00 



00 

«=*:+! 



^/2 S 



(n-fc)/2+!i/4 



(5.6) 



The lowest order surface coefficient coincides with the Dirichlet one, 



bi/2 = 6f/2 = -"V" 



(5.7) 



whereas the higher order ones take the following form in view of the known expressions 
for the coefficients of the nonanalytic and analytic parts of the action, 



bn/2 — 



■ r(fc-2j + i-|) ^^^^2,-1 

r(i)r(i-j) 
^r(i-fc + |) 1 



J > 1, 



^1-5/2 



r(^ + |) (t 



i~ky. 



>k + i. 



The first few coefficients bj are ultraviolet finite, 

J V r(fc)r(|-j) ^ - J- 



and 



(5.8) 
(5.9) 

(5.10) 
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However, for 2j > k they are divergent, 

but the complimentary coefficients bj have poles with residues that are exactly opposite in sign, 
{bj + bj)P°^'^ = 0. This is certainly a manifestation of the cancellation of infrared divergences between 
the analytic part of p.lOp related to bj and the non-analytic part related to bj . As a result the heat 
kernel stays well defined but acquires logarithmic terms in s because 

bj + bj s'-^/^ = ( (3j ln(W) + 7j ) , (5.12) 

3- - 2^-'^H XV^^ r(2j-2) ^2.7-1 ir. -|o^ 



The rest of the coefficients with half-integer numbers &j+i/2 are finite, 

m 

2T{2j + l-k) T{k) 



Thus finally the heat kernel trace in even integer (unregulated) dimensionality of a branc d ^ 2k 
takes the form 

Trif+^^ e^^f-^+i) = - f d'^+^X 

(47rs) — 



[k/2] 



oo I 

Y {(3, Hsm')+lj)s' \- (5.15) 



2j>k+l 

In odd dimensions, d = 2A;-|- 1, the poles of infrared nature are absent, so that both bj and 6„/2 are 
finite, and no logarithmic terms arise in the heat kernel. One can check that its trace has explicitly 
the form (|5.6p with d = 2fc -I- 1. Thus it has not only half- integer powers of the proper time, but also 
powers multiple of a quarter. 

Note that 61/2 for all finite m has the value (|5.7p characteristic of the Dirichlet problem. At 
the same time, for m — > 00 our problem reduces to the case of Neumann conditions, cf. Eg. p.lOp . 
corresponding to = |17j . However, this does not present any contradiction, because the 

asymptotic expansion in 1/M is obviously not homogeneous in m — > cx), as it involves growing positive 
powers of m. To analyze the limit of large m or e ^ 00, which includes the case of the massless DGP 
model, M — 0, we have to consider another representation of the effective action. This is discussed in 
the next section. 

6. Massless limit: effective potential in a toy DGP model 

As shown in Appendix B, the brane part of the effective action (j2.7p can be exactly represented in 
terms of the hypergeometric function. This representation has the form 



2 '""^ dTid/2)J "''sinf [att J '^^ 

'^Jd^-Y -^I'^^l-^'l--')' (6-1) 



2 

'aPV r 



2 

2\d/2 1-d 

4^ I ~2d^ 
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where F(a, b; c; u) is a hypergeometric function given by (jB.4p . and its argument u is defined in terms 
of the following two functions of e = m/2M, 



v± 



= ±y/l + e^ -e = -l/v^. (6.2) 



This representation exactly recovers the nonanalytic part of the 1/M-expansion (|3.10p . whereas the 
analytic part of the latter follows from the known hypergeometric series F{a,b;c;u) = 1 + 0{u) in 
powers of the argument u = 1 — w| ~ e. 

In the opposite limit of large e (corresponding to large m or small M) the brane action can be 
equivalently represented in a form useful for the expansion in u+ '--^ l/2e ^ 1. As shown in Appendix 
B by using the transformation formulae for F{a,b; c; u) it reads 

r2^ 



I Tv^d)^^ phr.no ^ j d-^X |-ir(-d/2)(l-Z;2)''/2 

TT (1 - «2 )d/2 



2r(l + I) sin(7rd) 



2) 0111(^/1 uy W^J. 



r - 



1 3 + d 
2' ~ 



3 



{d+l)F[l, ^;vl 



-F[l, ^: — vl 



(6.3) 



Using this representation one can consider the limit of exactly Neumann boundary conditions 
corresponding to m — > oo and a finite value of the mass M in the bulk. In this limit ii^ — s- and 
the last two terms of (|6.3I) vanish. Naively, the second term ~ e'^ is growing to infinity, but in the 
domain of ultraviolet convergence d < 0, so that with the appropriate order of taking the limits (first 
in TO — > oo and second in the dimensionality) it vanishes, too. Thus, only the first term remains, 

d 

- Tr^'') lnF^™'° = r(-d/2) / d'^x, (6.4) 



2 4 ' ' ' V 47r , „ 

which is obviously one-half of the contribution of the d-dimensional massive particle corresponding to 
the brane-to-brane mode propagating with the square-root operator jr'i="^"'= = y/ — □, 

i Tr^'^) In f'^™'° = i Tr''') ln(M2 _ □). (6.5) 

When added to the Dirichlet effective action (|2.4p it alters the sign of the brane (d-dimensional integral) 
term, which corresponds to the transition from the Dirichlet value of the surface coefhcient (|5.7p 
to the Neumann value, 

- (6.6) 

Another interesting limit of e — > 00 corresponds to the massless case of the DGP model with 
AI = 0. In this case the only nonvanishing term is contained in the second line of (|6.3p . and it yields 



1 Tr(^) InF^— = ^ ( ^ ) , j' . , / d'^x, (6.7) 



2 \4:TT J r(l + f) sinTTd 

because M/v+ — > 2eM — m. This result can be confirmed by a direct calculation of the effective 
potential of the brane mode, circumventing the operation of taking the limit M — s- in the general 
answer (16.31). 
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Indeed, the eflFective potential V^ff(m), 

i Tr^'') InF'^'""" = J d'^xV,s{m), (6.8) 

for the brane operator in the toy DGP model, iT'trano _ [—D + my^—\3)/m, can be written down in 
the form of a momentum-space integral similar to (|2.7p . 

oo 

Vctf{m.) = x^/,-r./rfN / dpp'^^^ \n(p^ +mp) 



(47r)'i/2r(|) 

u 

(Sr^FTi) /^"'"' [ln(t+l)+lnt + lnm2] . (6.9) 



Within the dimensional regularization only the first logarithmic term gives a nonvanishing contribution 
which we transform by using the proper-time representation of the logarithm and changing the order 
of integrations. 



^ ' V4^y r(|) 



oo 





1 / \d/2 1 



47r / r(| + 1) sinTTd' 



(6.10) 



After subtracting the ultraviolet divergence in the limit of the physical dimensionality, d N, 
this gives the renormalized effective potential of the usual Coleman- Weinberg structure. 



1 (-1)^ /m2\^^/2^^^m 

2 r(^ + l) V 4^/ 'V 



yM--.7^nr^(Tz) In-. (6.11) 



Here /x is the parameter reflecting the rcnormalization ambiguity, and the role of the field (argument 
of the potential) is played by the scale m of the brane term in the classical action (|1.7p - the analogue 
of the DGP scale. This result confirms the boundary effective action calculation of [S]. 



7. Conclusions 

In conclusion, we have derived the one-loop effective action in a simplified model of brane-induced 
gravity which gives rise to special boundary conditions involving second order tangential derivatives. 
The main peculiarity of this action is the presence of logarithmic ultraviolet divergences for both even 
and odd dimensionalities of the spacetime. This is different from analogous one-loop calculations in 
spacetimes without boundaries leading to divergences only for even spacetime dimensionalities. The 
action has a nontrivial analytic structure in the mass parameter — for generic M and m it is given 
by two representations in terms of hypergeometric functions (|6.1I) and (|6.3p and simplifies for the case 
of the massless field in the bulk, M = 0, to the form (|6.10p . After ultraviolet renormalization this 
gives rise to the familiar logarithmic Coleman- Weinberg effective potential ~ (p'^ \n{(p'^ / fi^) with the 
field (fi played by the parameter m in the boundary conditions (jl.lOp — the result used in [5] for the 
stabilization of the DGP cross-over scale in the Brans-Dicke modification of the DGP model [77] . 

We also derived the proper time expansion for the functional trace of the heat kernel subject to 
these generalized Neumann boundary conditions. This turns out to be nontrivial, because for even 
dimensionalities of the boundary it involves together with the well-known half-integer powers of the 
proper time s also the logarithmic terms ~ Ins, cf. Eq. (|5.15p . and for odd dimensionalities contains 
powers of s which are multiples of one quarter, see Eq. (|5.6p . Such peculiarities of the proper time 
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expansion are usually associated with the lack of strong ellipticity of the boundary value problem 
pi] when a naively positive elliptic operator acquires due to the presence of the boundary an infinite 
set of negative modes. These modes make the heat kernel operator unbounded and violate usual 
assumptions underlying its proper time expansion. But, as we show in Appendix C, strong ellipticity 
of our problem gets violated only for negative m in ()1.10|) — when the relevant classical action ()1.7p 
and the brane-to-brane operator p^'^^'^^^ Eq. (jl.lip . are both not positive definite. However, the heat 
kernel expansions become exotic also for m > with no violation of strong ellipticity, which implies 
deeper reasons of these peculiarities. 

To summarize, we conclude that the technique for quantum effects in brane models is more com- 
plicated than in systems without boundaries. Moreover, it does not reduce to a simple bookkeeping 
of surface terms in the heat kernel expansion of [T71 [THl [H] , and so on, because of the complicated 
square-root structure of the brane propagator (|l.lip mediating the effect of the generalized Neumann 
boundary conditions (jl.lOp . The proper time method that was fundamentally efficient in models with- 
out boundaries [14i il6j in our calculations above became a derivative of an alternative calculation. 
Namely, the surface terms in the heat kernel expansion were recovered from the 1/M-expansion of the 
action obtained by a different method of a Fourier decomposition. 

Nevertheless, the proper time method still does not loose its power and can be used in realistic brane 
models including gravity. In these models the effective action should be expanded in powers of the 
bulk spacetime curvature and the extrinsic curvature of the brane, starting with the approximation 
considered above. The momentum space decomposition is not very efficient for sake of such an 
expansion — the difficulty usually circumvented in background field formalism by the use of the 
Schwinger-DeWitt proper time method |141 1151 116j . Here we present without derivation (that will be 
given in forthcoming publications) such a representation for the Green's function of the DGP brane- 
to-brane operator and its functional determinant. They have a form of the weighted proper time 
integrals 

oo 

dse'°w{s), (7.1) 



OO 

Trln(-□ + rn^/^)=-Tr /^e^°i±|^, (7.2) 





2 



with the weight function u'(s) given in terms of the error function <f>(x) — dy exp(— y^) and 

having the following ultraviolet and infrared asymptotics 

The advantage of this representatiorQ is that it applies also to the case of the curved-space d'Alember- 
tian □, so that the generally covariant expansion of (|7.ip - (|7.2p in curvatures can be directly obtained 
by using a well-known Schwinger-DeWitt expansion for e'*'^. Thus the lowest order approximation 
for the exact brane-to-brane operator (jl.Sp in models with a curved bulk and curved branes can be 
considered by means of the manifestly covariant technique which can be systematically extended to 
higher orders. Combined with the method of fixing the background covariant gauge for diffeomorphism 
invariance in brane models, developed in [28], this will ultimately give the universal background field 
method of the Schwinger-DeWitt type in gravitational brane systems. 



^Note, in passing, that the interpretation of this weight and its asymptotics is very transparent. In the ultraviolet 
domain of small proper time ruy/s <C 1 (or big 2> m) it approximates the brane operator by — □, whereas in the 

infrared domain niy/s 3> 1 (or V— D <S m) it corresponds to its low-energy behavior m\/— CI. AH the results above 
could be obtained with the aid of this representation generalized to the case of a nonzero M. We did not use it, however, 
because this generalization has a complicated weight function. 



13 



A. Asymptotic expansion for integrals with a weak singularity 



The analytic and nonanalytic parts of the asymptotic expansion for the integral ()3.4|) (a so-called 
integral with a weak singularity) can be found by the method of [2Bj. For an integral of a slightly 
more general form this expansion reads 

no 

T(—a) n'. 

OO 

+ ^a„e", e^O. (A.l) 

n=0 

Here (p(t) is a function which is analytic a.t t — and has a Taylor series with coefRcients (y5^"''(0) = 
cP if / dt^ {Q) . The parameter (3 is positive in order to guarantee the convergence of the integral at i = 0. 
The first sum gives a nonanalytic part of the expansion determined entirely by the derivatives of the 
function ip, (/^("^(O), whereas the second sum determines the analytic part with coefficients involving 
a (nonlocal) dependence of the function ip{t) at all t. These coefRcients are given by the following 
expression: 



rpa+f3+m-n (pi™) (Q^ 



m=Q 



where (ps{t) is a piecewise smooth function obtained from ip{t) by subtracting its first few terms of 
the Taylor expansion at < = on a finite segment of the t-axes, < t < T, 



^s{t) = ^{t)- ^ ^-iit™, t<T, 



ml 

m— 

ps{t)^^{t), t>T. (A.3) 

Here the number of subtracted terms is given by [n — a — (3] — the integer part of rt — a — /3, and T is 
arbitrary positive. The value of the latter is immaterial, because it is easy to check that 9a„/9T = 0. 
These subtractions are necessary to guarantee the convergence of the integrals in (|A.2p at t = 0. For 
the first few a„ these subtractions are absent, 

oo 

a(a — l)...(a — n + I) f „ , ,. ^ 

an = — / div? i n<a + f3, A.4 

n! J 



while for n > a + /? their effect can be explicitly circumvented by multiple integrations by parts 
in (|A.2|) . After {Nn + 1) integrations by parts, all non-integral terms of (|A.2p cancel out and the 
expansion coefficients take the following alternative form, 

_ a{a ~ l)...{a - n + 1) T{n ~ a ~ (3 - Nn + I) 
^" ~ n\ T{n-a- [3+ 1) 

OO 

1 



a-(3-Nn 



^(^"+i)(t)r+^-"+^", n>a + (3, (A.5) 



where 

N„ = [n-a-(3]. (A.6) 
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Finally, when a + /3 is a positive integer N the two sums of (|A.1[) formally become analytic, but their 
coefficients develop pole singularities in a + /3 — iV — > 0. These singularities come from these two sums 
with opposite signs and cancel. The finite remnant of this cancellation is a term logarithmic in e. 
If we apply now this asymptotic expansion to the case of our integral (|3.3[) - (|3.4[) with 



ipit)^iVlT¥ + tf^', a = (A.7) 



then in view of 



= (-1)"— TTTTT^T ' (A- 



r(-a) ^ ' sin(7rd/2) 

the nonanalytic part is explicitly expressed in terms of the function ip{—s) with the flipped sign of its 
argument. 



sin(7rd/2) n! sin(7r(i/2) 



The coefficients of the analytic part can be explicitly calculated by taking the integrals (|A.4p in the 
domain n < d/2 — 1 and extended beyond this domain by analytic continuation (which is equivalent 
to using ()A.5|) ). The result is 

Therefore for odd n ~ 2j + 1 they vanish due to the unregulated (by the dimensionality d) infinity in 
the denominator, whereas for even n = 2j they are given by p.9p . and the final form of the expansion 
for / is given by p.7p . 

B. Hypergeometric function representation 



The integral ()2.7p can be rewritten in terms of the integration variable x — ^/p^ + A'P /M as 

oo 

j dpp'^^^ In [jp- + my/NP 



r2 4_p2 



= M^/.„,.>.l,^/»,.(.^ + 2..-l). (B.1) 

1 

By factorizing the argument of the logarithm and integrating the result by parts we convert it to the 
sum of two terms. 



j dxx {x^ - 1^ (^^2 - l) 

1 

oo 



oo 

1 J2 J dx {x' - If/^x - v)-\ (B.2) 



d 

v = v± 
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where v± are the roots (|6.2p of the quadratic polynomial x"^ + lex — 1 = [x — Vj^)[x — vJ). 
With the change of the integration variables x = we have 

oo 1 

dx {x" - Xfl'' {x - v)-' = 1 Jdt ^-^-''/'(l - (1 - tv^r' 



1 

+ 1 J dt t-l/2-d/2 (1 „ ^)<i/2 (1 _ ^^2yl^ (3 3) 



SO that finally in terms of the hypergeometric function 

F{a,b;c;u)^—^^j^^— f dtt''-'^ {1 ~ ty-''-^ {1 ~ tu)-" (B.4) 



T{b)T{c-b)Jo 
the basic integral reads as 



dx (x^ - if/^ (x - v)-^ = fm/'^ ^('' 



V r(i + |)r (i^) i-d. 3 



F l, ^— . (B.5) 



2 r(f) 

We will need its ^-expansion. At e — > the parameter v — v± ±1, so we have to transform 
the hypergeometric functions to the series in 1 — v^. Because of the known relation F(a,b; a;u) = 
F[b, a; a; u) = (1 — u)^'' we have 

F[l,-^;l;v')={l-vy/' (B.6) 
and in view of Eq. 9.131.2 of f26l, 



F 1 



V ' 2 ' 2' J r(i)r(i + f) V ' 2 ' 2' 



^ ' r(i)r(i^) V2' 2' 2' /' ^ ' 



where again f( i, 1 + |; 1 + |; 1 - f^) = [1 - (1 - v^)]-^/^ = \v\-^. Therefore 



F I, ——■ -■ V 



2 ' 2' " y ' d T{^) 



2 

Substituting (|B.6p and (|B.8|) into (|B.5|) we have 

OO 

dx {X^ - „ ^ ^^(^;) (1 _ y2^d/2 

sm(7ra/2) 

where a step function arose as the result of summation of two terms, 
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The first term in (|B.9[) exists only for positive v and is nonanalytic at v ^ 1, whereas the second 
term is analytic. Obviously for negative v this integral is an analytic function because the argument 
of the logarithm nowhere tends to zero in the integration domain. This explains the absence of the 
first term for v < 0. 

Substituting (|R9|) to (|RT|) - (|R21) we finally get 

OO 



V — V± 

Bearing in mind that (1 — f^)*^/^ ~ {2sv^y^'^'^ ~ {2eY^^(f{—e) we finally have the representation useful 
for small 

oo 




which gives rise to the representation (|6.ip . 

To consider the limit of e — > cxd we need another representation, because in this limit v+ ~^ l/2e 
and V- ^ — 2e — cx3, so that the contribution of u = in (IB.2P should be expandable in 1/v. This 
can be achieved by the transformation formula 9.132.2 of [26] which in our case reads as 

/ l-d 3 ,\ 1 1 / 13 + dl\ 

+ 2 r(l + |) (-i;2)i/2 ■ ^"-^^^ 

Using this in the representation (|1.6p of the v — v- term of (|B.2[) and taking into account that 
l/u_ = — f+ we finally arrive at the equation underlying the representation (j6.3p : 



^ dp/"^ In + my'AP +p^^ 



2d sin(7r(i/2) d sin(7r(i) 

1-d 3 



4v^ 



(d+l)F(l, 

-F(l,l;^;4 
2' 2 ' + 



(B.14) 



C. Strong ellipticity problem 



The strong ellipticity problem for the operator p.Sp with generalized Neumann boundary condition 
(|1.10P consists in the existence of an infinite set of normalizable eigenmodes with a spectrum which 
is unbounded from below [24| . This implies that the relevant heat kernel is an unbounded operator 
which cannot be rendered bounded by the elimination of a finite number of states from its functional 
space, and therefore it has an unusual proper time asymptotics different from p.l2p . This situation 
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occurs when the parameter m in the operator (|1.9p of the boundary condition is negativ^l. Here 
we show that these negative modes correspond to the negative modes of the brane-to-brane operator 
(jl.lip . locahzed near the brane/boundary and responsible for brane-to-brane propagation. 

Indeed, in this case there is a set of eigenmodes locahzed near the boundary y = of the form 

<^p{x,y) = e'^^-^-y (C.l) 

in which Ap is given on account of the boundary condition by the expression 

2 

Ap = -— > 0, TO< 0. (C.2) 
m 

Negative m guarantees the normahzabiUty of these eigenmodes which, therefore, cannot be excluded 
from the functional space of the operator. Their eigenvalues Ap, 

{M^ -n^a+i))%{X)^Kj,%{X), (C.3) 

Ap=(M2+p2_i^), (C.4) 

are negative for sufficiently high Fourier momenta 

2 

Ap<0, p2 > ^ (l + Vl + 4M2/m2), (C.5) 

and tend to — oo for —> oo. But the momentum space domain where they are negative exactly 
coincides with the domain in which the brane operator is negative definite for m < 0, 

F""^--- ^pix) - + v/M2+p2^ ippix), <^p(x) = e'P-. (C.6) 

Thus, the lack of strong ellipticity of the generalized Neumann boundary value problem is in fact 
the lack of positivity of the action (jl.7p with m < 0, from which this problem originates by the 
variational procedure. 
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